function [dat]=r_eke(year)
% ---- Read EKE parameter
% (from DUACS data)

dat = struct('dummy',1);

dir_duacs='/calcul/pc-charria/OBSERVATIONS/SATELLITE/SLA/';

D=dir([dir_duacs,'*.xyz']);

n=0;
for il=1:length(D)
%for il=271:274
    if (strfind(D(il).name,year))
        n=n+1;
        dd=load(['/calcul/pc-charria/OBSERVATIONS/SATELLITE/SLA/',D(il).name]);
        sla(n,:,:)=reshape(dd(:,3),519,470)./100;
    end
end


dat.londata=unique(dd(:,1));
dat.londata(find(dat.londata > 180))=dat.londata(find(dat.londata > 180))-360;
dat.londata=sort(dat.londata);
dat.latdata=unique(dd(:,2));

        

% [dat2]=load_gadget([1997 (2000+0.9999999999)],[256 52 -21 81.2],'ssha_ddmr_v10.nc');
% 
% dat.londata=dat2.lon;
% dat.londata(find(dat.londata > 180))=dat.londata(find(dat.londata > 180))-360;
% dat.latdata=dat2.lat;
% sla=squeeze(dat2.ssha);
sla=permute(sla,[3 2 1]);
% 

% ... Warning: Remove the equatorial band !!
eke_tot(1:size(sla,1),1:size(sla,2))=0.;
for it=1:size(sla,3)
a=squeeze(sla(:,:,it));
%a=nanmean(sla,3);

%if (DUACS == 1)
% 1D case
a(:,find((dat.latdata > -2) & (dat.latdata < 2)))=NaN;
%else
%    a(find((lat > -2) & (lat < 2)))=NaN;
%end
res_eke=eke(dat.londata,dat.latdata,a.*100.);

eke_tot=res_eke+eke_tot;
end


dat.data=eke_tot./size(sla,3);

% ---- Clean the Structure ...
dat=rmfield(dat,'dummy');    
